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TECHNICAL NOTE 2873 


THE EFFECT OF LONGITUDINAL STIFFENERS LOCATED ON ONE SIDE 
OF A PLATE ON THE COMPRESSIVE BUCKLING STRESS 
OF THE PLATE-STIFFENER COMBINATION 


By Paul Seide 


SUMMARY 


The problem of buckling under uniform compression of flat, simply 
supported, rectangular plates with equally spaced longitudinal stiffeners 
on one side of the plate is investigated. For the case of a plate with 
one, two, or infinitely many stiffeners, the analysis yields expressions 
for the effective moment of inertia of the stiffeners that can be used 
in (© nią with the buckling charts previously presented in NACA 
TN 1025 o 


INTRODUCTION 


The buckling of stiffened plates is & subject which has received much 
attention in the literature on aircraft structures. With few exceptions, 
the solutions presented are idealized and are valid only for stiffened 
plates for which the center of gravity of each stiffener cross section 
lies in the middle surface of the plate or for which the stiffeners are 
hypothetically connected to the plate in such a manner that sliding of 
the stiffener along the plate surface is permitted. Although these solu- 
tions clearly establish the relationship between the buckling stress of 
the plate stiffener combination and the flexural stiffness of the stiffener, 
they do not determine the effective flexural stiffness provided by stif- 
feners in the usual aircraft application, that is, when they are riveted 
to one side of the plate. 


Timoshenko suggested (ref. 1) that the stiffness of a one-sided 
stiffener might be taken into account by replacing the moment of inertia 
of the stiffener about its center of gravity by an effective moment of 
inertia and using this value in the solutions valid for stiffeners with 
their centers of gravity located in the plate middle surface. In several 
illustrative examples he took this effective moment of inertia as the 
moment of inertia of the stiffener cross section about the plane of con- 
tact with the plate. That this method of correction is arbitrary and 
not generally applicable to all plate stiffener proportions, however, can 
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be readily seen. If the stiffeners are very large compared to the plate, 
the attached plate can have very little effect on the bending of each 
stiffener about an axis passing through the center of gravity of the 
stiffener cross section. On the other hand, stiffeners that are very 
small compared to the plate would appear to be forced to bend about the 
plate middle surface. The effective moment of inertia of a given stiff- 
ener can therefore vary between the moment of inertia taken about the 
center of gravity of the stiffener cross section and the moment of inertia 
taken about the plate middle surface depending upon the proportions of 
the attached plate. For Z-stiffeners of the proportions encountered in 
aircraft construction, the ratio of the moments of inertia based on these 
two limiting positions for the assumed neutral axis of bending is greater 
than 2.5 to 1. Because the effective moment of inertia chosen to repre- 
sent a given stiffener strongly influences the calculated buckling stress 
of the plate stiffener combination, it is therefore desirable that any 
arbitrariness in its calculation be eliminated. 


Investigations of plates with stiffeners on one side have been made 
and verify the dependence of the effective moment of inertia of the stiff- 
ener on the relative dimensions of the plate and stiffener. These solu- 
tions, however, are generally limited to special cases. The bending of 
an infinitely wide plate with a longitudinal stiffener on one side is 
investigated in references 2 to 4. In reference 2, Timoshenko and 
Goodier neglect the flexural stiffness of the plate and consider loads 
applied only to the stiffener; an extension is made in reference 3 by 
Smith, Heebink, and Norris to include the plate flexural stiffness; both 
plate flexural stiffness and loads applied to the plate are considered 
by Odqvist in reference 4. Shear buckling of a finite-width plate with 
a centrally located longitudinal stiffener on one side is investigated 
by Chwalla and Novak in reference 5. Good accuracy would be expected in 
their analysis only when the plate is relatively wide. Cox and Riddell 
(ref. 6) give an expression based on an effective-width concept for the 
effective stiffener stiffness for compressive buckling of a finite width 
plate with a centrally located longitudinal stiffener on one side. The 
stability of plates with closely spaced longitudinal and transverse stiff- 
eners is investigated approximately by Freyer in reference 7. A very 
accurate set of differential equations for the same problem is developed 
by Pfluger in reference 8. These differential equations contain terms 
which reflect the consideration of changes of geometry of the plate 
beyond those considered in the elementary theory of bending and buckling 
of flat plates. The contribution of these terms is shown to be insignif- 
icant, however, in the practical range of stiffened-plate construction. 


The references cited indicate that some studies have been made of 
the buckling behavior of plates with one stiffener or with many closely 
spaced stiffeners. No material appeared to be available for the cases 
that fall between these extremes. In the present paper, a unified analy- 
sis is made of the buckling behavior under uniform compression of a flat, 
rectangular, simply supported plate with one, two, three, or infinitely 
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many longitudinal stiffeners located on one side of the plate. The 
essential feature of the results of this analysis is that the effective 
moment of inertia of one-sided stiffeners on plates of arbitrary dimen- 
sions can be written in terms of a simple correction to the moment of 
inertia of the stiffeners about their own centers of gravity. This 
effective moment of inertia can then be used in conjunction with the 
buckling charts of the type previously presented in reference 9. The 
analysis leading to this result is given in the appendixes. 

The material presented herein was submitted to the University of 
Virginia in partial fulfillment of the requirements for a masters degree. 


SYMBOIS 

Ag stiffener area 

a plate length 

an Fourier coefficients in bending deflection function 

b plate width, Nd 

Et? 
D plate flexural stiffness, 5 ‘ 
12(1 - | 

d bay width 

E Young's modulus for plate and stiffeners 

Fy shear force per unit length applied at ith contact line 
for plates with one, two, or three stiffeners 

Po shear force per unit length applied at any contact line 
for infinitely many stiffeners 

Fin coefficients of shear forces for plates with one, two, 
or three stiffeners 

Fn coefficients of shear forces for plates with infinitely 
many stiffeners 

J,8 integers 

i integer denoting contact line or bay 

k number denoting location of contact line with respect to 


plate center line 


zo aa 
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moment of inertia of stiffener cross section, taken about 
axis parallel to plane of plate and passing through 
center of gravity of stiffener cross section. 


effective moment of inertia of stiffener 


integer denoting number of buckle half-waves in longi- 
tudinal direction 


number of bays 


additional forces per unit width in plane of plate 
middle surface in ith bay 


integers denoting mode of buckling 


additional load in ith stiffener 


integers denoting number of buckle half-waves in trans- 
verse direction 


plate thickness 
total potential energy of stiffened plate 
middle-surface displacement 


additional displacements in plane of plate middle surf&ce 
in x- and y-directions in ith bay, respectively 


additional displacements in x-direction at center of 
gravity of cross section of ith stiffener 


strain energy 

potential energy 

bending deflection, normal to plane of plate 
distances along coordinate axes (see fig. 1) 


modal coefficient appearing in expression for effective 
moment of inertia of stiffeners 


distance normal to plane of plate between plate middle 
surface and center of gravity of stiffener cross section 


additional direct and shear strains in plane of plate 
middle surf&ce in ith bay 
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Eg | additional direct strains at center of gravity of cross 
i section of ith stiffener 
1137127713; modal coefficients appearing in criterion for buckling 
symmetrical about longitudinal center line of plate 
742715716 with three stiffeners 
TI Poisson's ratio for plate material 


ài yt buckling-stress parameter 
xD 
a/a buckle aspect ratio 
a/a aspect ratio of each bay 
-— flexural-stiffness ratio (called y in ref. 9) 
ett effective-flexural-stiffness ratio 


area ratio 


— — a my 9— € 
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Subscripts: 

N number of bays 
P plate 

S stiffeners 


RESULTS AND DISCUSSION 
Statement of the Problem 


The problem that is investigated in the present paper is the 
determination of the buckling stress of a flat, simply supported, 
rectangular plate with one, two, three, or an infinite number of longi- 
tudinal stiffeners on one side. The stiffeners are equally spaced, are 
each of the same cross section, and have a Young's modulus equal to that 
of the plate. (See fig. 1.) The plate is loaded by uniform compressive 
forces per unit width o,t on its transverse edges and the end load on 


each stiffener is a compressive force O Ag. 


Assumptions and Limitations of Analysis 
The assumptions and limitations of the analysis are as follows: 


. (1) Each stiffener is assumed to be continuously and rigidly con- 
nected to the plate along a straight line (the contact line) in the 
plate middle surface. This condition may be visualized by supposing 
the stiffener to be extended to the contact line by a rod of infini- 
tesimal area. (See fig. 2.) The effects of rivet spacing and of rivet 
flexibility are thus excluded from the analysis. 


(2) Each stiffener is assumed to have zero flexural stiffness for 
bending parallel to the plane of the plate and zero torsional stiffness. 


(3) Small-deflection theory for bending of elastic plates and 
elementary beam theory for bending of the stiffeners are used. The 
analysis is thus limited to a general instability of the plate-stiffener 
combination and excludes such effects as local buckling of stiffeners, 
cross-sectional distortions, and shear deformations of the stiffeners, 
all of which decrease the effective moment of inertia of the stiffener. 


(4) An artifical boundary condition is introduced into the analysis; 
that is, points in the transverse edges of the plate are taken to be free 
to move transversely before buckling occurs but are completely restrained 
from further movement after buckling. Any other boundary conditions would 
needlessly complicate the analysis. The longitudinal edges of the plate 
are taken to be free from middle-surface forces both before and after 
buckling. 


SSS rE E — ———— —— m E — 9 — m mać 
Ss 
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Summary of Analysis 


The stability criterions for compressive buckling of longitudinally 
stiffened plates with deflection of the stiffeners, derived in appendix A, 
with the exception of the criterion for buckling symmetrical about the 
longitudinal center line of a plate with three stiffeners, can be 
expressed as follows: 


o A fafa a ot ^ 
Els, 1 EAgZg x2 \ m x Pi z (2/4) Ag 
dD 1+ Zya 28 aD sin 04 sinh 0o xD NM dt 
os ai 02 
cos = cos 01 cos T- cosh 09 
(1) 
where 
O. = 
1 
(2) 
05 = 


(m = 1, 2, 35 . . 00) 


and. Za is a modal coefficient. This modal coefficient ís a function 


of the buckle aspect ratio = the number of bays N, and the mode of 
buckling indicated by q as in the following table: 


Number of Mode of buckling 
stiffeners, about longitudinal 
center line of plate 


Symmetrical 
Symmetrical 


Antisymmetrical 
Antisymmetrical 
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d 
The variation of ZNq with the buckle aspect ratio xi is shown in 
table I and in figure 3. Poisson's ratio | was taken as 1/3 in the 
calculation of these values. 


The criterion for buckling symmetrical about the longitudinal 
center line of a plate with three stiffeners is given by 


1 2B -2 o 2 
O SA S Ag _ 
aD 2 aD ne m at 
A D 
e Zl: 59) - (ni E) 
2 dt 3 dt dt 
3 la^s.t 
n a/a x Ag 
ZA WET L Le => -_2 
ARN Es, —  — — "*"5dt ——  — EAszs | 
sin 04 sinh 05 dD Ag Ag Ag 2 dD 
a. e So, 

1 ] 2 2 dt 3 dt dt 

V2 cos 01 Y2 - cosh 5 
Lu 
2. 2 + EBI | e R 
d sx (a/a) åg ye Nm xD 
95 01 
Ye, cosh 05 ve + cos 04 
2 
Ag 2D 
16 dt EAgZg 
6 — | =0 (3) 


2 dD 
A 
(1217 SI 13 ÓN Es) l (m = l, e, E 6 «= co) 


22 
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d 
where 74 to ng are modal coefficients and are functions of E. 


The variation of Th to "le with zi is given in table II and is 
shown graphically in figure +. 


When buckling symmetrical about the longitudinal center line of a 
plate with one, two, or infinitely many longitudinal stiffeners on one 
side is investigated, the buckling charts of reference 9 may be used. 


In place of the flexural-stiffness ratio = these charts are entered 


2 


with values of the effective-flexural-stiffness ratio 


_ 2 
Eleff Elg T EAgZg (4) 
aD db e" Ag dD 
“Na gt 


where EI, ep/dD is calculated with the proper value of the modal coef- 
ficient Zug obtained from table I or from figure 3 of the present 


paper. Thus, if the buckling-stress parameter for a certain plate- 
stiffener combination is desired, the chart of reference 9 for the 
appropriate area ratio and number of stiffeners is entered with 
ETofp/dD and the given value of the plate-bay aspect ratio to obtain 


the buckling-stress parameter. In this procedure, the value of m used 
in reading Na from figure 3 of the present paper must correspond to 


the value of m in the region of the charts of reference 9 from which 
the buckling-stress parameter is read. In the charts of reference 9, the 
value of m in the region to the left of the first dashed-line curve 

is l and the value of m in the region to the right of that curve is 2; 
m would increase in value by 1 in the region to the right of each 
succeeding dashed-line curve. 


When buckling symmetrical about the longitudinal center line of a 
plate with three stiffeners on one side is investigated, equation (3) 
must be solved by trial and error for the value of the buckling-stress 
parameter corresponding to a given value of m. Different values of m 
are tried until a minimum value of the buckling-stress parameter is 
obtained. 


Investigation of the antisymmetrical modes of buckling for plates 
with two or three stiffeners involves a trial-and-error solution of 
equation (1) for the buckling-stress parameter for given physical 
dimensions of plate and stiffeners. 
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For each case the criterion for buckling with nodes at each 


stiffener 
c 2 
d ot d. 
x - (a E) (5) 


nod a/a m 


should be checked since this criterion defines the highest buckling- 
stress parameter obtainable in conjunction with an assumption of zero 
torsional stiffness for the stiffeners. 


CONCLUDING REMARKS 


The buckling under uniform compression of flat, simply supported, 
rectangular plates with equally spaced longitudinal stiffeners each of 
the same cross section and a Young's modulus equal to that of the plate 
is investigated. The effect of one-sided placement of stiffeners on the 
plate is incorporated in a simple expression for the effective moment of 
inertia of the stiffeners. For a plate with one, two, or infinitely many 
stiffeners, the effective moment of inertia can be used in conjunction 
with buckling charts of the type previously presented in NACA TN 1825. 
For a plate with three stiffeners, a new stability criterion is derived. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., June 16, 1952. 
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APPENDIX A 


STABILITY ANALYSIS OF PLATES WITH ONE-SIDED STIFFENERS 


Determination of Stress Distribution in the Buckled Plate 


Consider the bending of a plate with stiffeners attached so that 
sliding of the plate and stiffeners is permitted.  Middle-surface dis- 
placements u in the plate along each contact line exist and in each 


stiffener there exist contact-line displacements ug + Zg Ss, where ug 
x 


is the displacement at the center of gravity of the stiffener cross sec- 
ow 
Ox 
sections of the stiffener being assumed to remain plane and perpen- 
dicular to the deformed axis of the stiffener. In general, u and 
ug + Zg QU 


tion and Zg represents the displacement due to bending, the cross 


are unequal and there would exist relative displacements of 


corresponding points of the plate and stiffeners. (In the present 


problem, for example, if the stiffeners were permitted to slide along 
«the sheet, u and ug would be equal but relative displacements due 


to bending, of magnitude Za A would exist.) Since in an actual plate 
x 


these relative displacements are not permitted, a system of shear forces, 
which produce additional deformations to cancel out the relative dis- 
placements, are induced in the structure. The system of shear forces 
consists of forces applied to the plate middle surface at each contact 
line and equal but opposite forces applied to the corresponding stiffener 
at each contact line. The effect of these forces on the buckling stress 
of a stiffened plate loaded in edge compression is investigated in the 
present paper. 


Prior to buckling, the stiffened plate considered herein is uni- 
formly compressed, the longitudinal direct stress in both the plate and 
stiffener being equal to -0,; all other stresses in the plate are equal 


to zero. When the stiffened plate buckles, however, self-equilibrating 
shear forces are induced at each contact line within the structure, as 
demonstrated previously, and the stress distribution within the plate 
changes although the applied edge stresses remain constant. The addi- 
tional stresses in the structure must be investigated before the buckling 
stress of the stiffened plate can be determined. 


- - — ar — — ey e A PP Bn — Hip M nl n ee —À C — 
-— m nr -1-— —.- s A nn — rs —— ee — ——— M mł nn — A Pn 
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Within each bay of the plate the additional middle-surface forces 
are continuous and in equilibrium, that is, for the ith bay 


Ney My 


Ox E dy 


(A1) 


L 
P gc (o, H s) (A2) 


du, 
*Xi Ox 

EI 
TS ma 


Xyi — dy Ox 


From the boundary conditions given in the section entitled 
"Assumptions and Limitations of Analysis," normal and tangential middle- 
surface forces must vanish along the longitudinal edges of the plate and 


RÁ, - — — n_n —À ee M 
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normal middle-surface forces and tangential displacements must vanish 
along the transverse edges of the plate; that is, at y = Nd 


2? 
Ny - Rr, = 0 (Alta) 
at y 2-2 
Ny, = Nay, = 0 | (Ab) 
and, at x=0 and x=a, 
Ny, 5V1=0 (1=1, 2, .. . N) (Alc) 


In addition, the following conditions must be satisfied along each con- 
tact line: 


Nu -N = F 
XJi XJ 41 i 
=N 
A Viel 
(A5) 
Ug 7 444 (1=1,2, . .. (N - 1)) 
Vi? Yi 


where F; is the shear force per unit length induced at the ith contact 


line after buckling. These conditions state that at each contact line 
the difference of the shears in adjacent bays is equal to the applied 
shear load, that middle-surface forces normal to the longitudinal edges 
of adjacent bays are continuous, and that middle-surface displacements 
are continuous. 


The additional forces and displacements along the centroidal axis 
of the stiffeners due to the shear forces induced after buckling are 
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found fram the equations of equilibrium of shear force and internal 
load to be 


dP; - 


E Fi (1ETL,2, s + « (No 1)) (A6) 


and the load-strain relationships are found to be 


P 
sa, = 35 (1=1,2, . . . (N- 1)) (AT) 
where 
dug 
ss, = (48) 


with the conditions that each load Pi must vanish at the ends of the 


corresponding stiffener. 


The appropriate expressions for the shear force at each contact 
line are now assumed to be 


Fį = Fig cos = [r1 2 w e (N21) (A9) 


where m is an integer. (It will be seen subsequently that this result 
is consistent with the assumption of sinusoidal buckling with m half- 
waves in the longitudinal direction and that this type of buckling 
actually occurs.) The forces and displacements in the plate and in the 
stiffeners due to a single shear force of the form of equation (A9), 
located anywhere in the plate, are derived in appendix B. The middle- 
surface forces and displacements at any point in the plate due to all 
the shear forces are obtained by superposing the individual forces and 
displacements at the point due to each of the shear forces. 


Thus far in the discussion the shear forces induced by contact- 
line restraint have been considered arbitrary. They can be related, 
however, to the bending deflections of the plate middle surface by 
conditions of continuity of plate and stiffener at each contact line. 
When buckling occurs, the addition displacements at each contact line 


in the plate due to the induced shear forces are ZĘ Ma and the 
U 
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contact-line displacements in the stiffeners are us, + sl SH x) Ma 
y=(1- 


In order for points in the plate and stiffeners to remain together, these 
two displacements must be equal at every point in the contact line. Then 


Ua +2 Ów = A10 
Si I Ma (alh. Pa (ALO) 


(1=1,2, . . . (N - 1)) 


An appropriate expression for the bending deflections of the stiffened 
plate is 


v= sin BED a oin Ñ SONS (A11) 


When equation (A11) and the displacements in the plate and stiffeners 
in terms of the coefficients Fim of the shear forces are substituted 


into equation (A10), & sufficient number of equations is obtained to 
give the coefficients F4, in terms of the coefficients a, of the 


deflection function. In general, these equations are of the form 


-1 
e (913 + Tijm)" jm = sal DE a, sin T (A12) 


j-l 


where O1 4 is the Kronecker delta (915 =] if i= j; 011 =0 if 
1 f 5) and. jm 7 jim and is given by the results of appendix B. 


Determination of Stability Criterions 


With the distribution of middle-surface forces and displacements 
in the buckled plate known in terms of the Fourier coefficients of the 
deflection function, the principle of minimum potential energy may be 
used to find the compressive buckling stress of the plate. The deflec- 
tion function given by equation (All) and the forces, strains, and 
displacements in terms of the coefficients of equation (All) are substi- 
tuted into the potential-energy expression for the structure. The 
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potential energy is then minimized with respect to each of the coef- 
ficients a. A set of simultaneous equations is obtained from which 


stability criterions can be determined. 


The potential energy of an elastic body is defined as the difference 
between the strain energy of the body and the potential energy of applied 
edge loads. At the onset of instability the structure has potential 
energy due to the applied end loads which remains constant while the 
plate bends and need not be taken into consideration. After buckling the 
plate acquires additional strain energy which can be expressed as the 
sum of the strain energy of bending of the plate and the strain energy 
of deformations in the plane of the plate middle surface (ref. 10) as 


a „Nd 2 2 e 2 2 
81 Is (Zw) «20 - w)|( 2) - 99 X|.ax dy - 
0 -Na Ox oy Ixe dy? 


a p(1-Da l N_ papri-Da 
EL iron die MAE HS EL Boda (ay, i 


N, e +N Y dx dy A13) 
Yi Ji xy, xy) 


The first term on the right-hand side of equation (A13) represents the 
strain energy of bending of the plate. The second term represents the 
strain energy arising from displacement of the middle-surface compressive 
forces, which remain constant as bending takes place, by additional 
stretching of the plate middle surface after buckling. The third term 
represents the strain energy due to stretching of the plate middle sur- 
face by the middle-surface forces that arise after buckling. 


The additional potential energy of the compressive forces applied 
along the transverse edges of the plate is 


- — loa u He - = 


AN a nk Be ml m B 


na ee ee —Á—À——— 


The first term on the right-hand side of equation (A14) is the work done by the edge loads in 
moving through displacements due to the shear forces induced after buckling. The second term 
represents work done by the edge loads in moving through displacements due to bending. 
of this term is consistent with the assumption in small-deflection bending theory that the plate 


middle surface undergoes no stretching due to bending. 


Similarly, the additional strain energy of the stiffeners is 


and the additional potential energy of the end loads on the stiffeners is 


2 


z N-1 > a 
| Wg S3 aia). œ Ms) o + "s =p J, CA Da dx 


The total potential energy of the stiffened plate is then 
ari EET AA Ki. da 
2 Jo J_Na dy ax? ay? i=1V0 (1-H-a)a dE 
i a (1- Da S a 2 
Id J gis Ct * Ty, ty, 7) 9 el [EP 5 
-1 pa -l pa i-a 
x Pigg, Y E GyÁgeg, dx + E. Se (n), Ñ (44) 29 da 
EST E, 092,] EO (Bgl S 


2 ME M 5 
dx + > | Pie SE € dx 
i*s xAgeg 
2 tio 1 isl V 0 i 


26 


tLoe MI VOYN 


The form 
(415) 
' (A16) 

5 


(A17) 


The potential-energy expression (eq. (A17)) can be considerably simplified by transforming * 
some of the terms in the equation. Consider, for instance, the term representing the work done a 
after buckling by the compressive forces applied y the D. in moving through displacements 


due to the induced shear forces - that is, S e O, , [e - (Ui) BL This term 
ial Yep xaQ 


"3 u 
can also be expressed as i f. 2) e.t i dx dy which is the equivalent of 
z Ñ x 

i=lJo V(1-5 -1)a ò 


(1-Da 
EIL = NE" dx dy and cancels the term representing the additional strain energy 
1=1 


in the s ius fram displacement of the constant middle-surface compressive forces after 


- - a 
buckling. Similarly, the term SH GzAg (4s ) = ("s ) cancels - S |, NS dx 
1 Xn"B5 i x-O| 1i=1Y0 t 


1=1 


The term representing the strain energy of stretching of the plate middle surface by the 
middle-surface forces that arise after buckling can be transformed to a much simpler expression 


by integrating the term by parts to obtain 
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The first, second, fourth, and fifth integrals on the right-hand side 
of equation (A18) vanish by virtue of the conditions at the plate edges 
(eqs. (A4)), the equations of continuity along each contact line 

(eqs. (A5)), and the equations of equilibrium of middle-surface forces 
within each bay (eqs. (A1)). The third integral, however, has a value 
and can be rewritten as 


a 
l 
i I (A m) yg 15 (A19) 
2 


But Fxyw is zero at y E and Nxy, is zero at y =-% so that 


the first and third integrals on the right-hand side of equation (A19) 
vanish. Furthermore, when the conditions of equations (A5) are con- 
sidered, the second integral can be rewritten as 


N-l pa 
> > y POD a ya dx (A20) 


which is the transformed additional strain energy of deformation of the 
plate middle surface. This result could also have been obtained by use 
of the law of conservation of energy; that is, internal energy equals 
external work. 


In a similar manner, the following transformation can be obtained: 


Nel pa Nel pa 
a SE 
2 J, Fig, dx 77g J Fius; E ag 
121 Jo 1=1JO 
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Equations (A20) and (A21) may be combined to yield 


r F; RZE: S S (A22) 


which, from the relationship given in equation (A10), can be written as 


i 
2 


ps 
I 
fF 


De S xm sg 


LH; 


The simplified potential-energy expression is then 


2 


o f SE A 2. 
U-2 2(1 - mde Co dx dy 4 
5 t J xa + 2( H) (e i a y 


Q/ 
no 


2 
N-1 ms [* (= ) 1 N-1 E dw i 
24 = Jo Nox reb ma mat t=1/0 CE =(1 Ma = 
t a pa” 2 N-1 2 
T i (=) dx dy - e zB , (8x) dx (A24) 


Substitution of the bending deflections w given by equation (A11) 
and the shear forces F; given by equation (A9) into equation (A24) 


yields 
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Ga” eo 2 [a/a e A Gztfa/af| o 
a U = > L AE Ss ol est an + 
nx ND n=1 yo\ mM xp NM 
2 2 N-1 2 
EL d c.t A = 
ee) 5 | (5 nai. 
N|dD ED e db | < (> ay sin » w 
EAcza’ N-1 — p co 
N aD 3 2. n | 
i=1 /m = n= N 
aym Z 
(=) os 


The third term on the right-hand side of equation (A25) is the contri- 


bution to the potential energy 


plate and stiffeners by the shear forces induced after buckling. 


of the structure of stretching of the 
If ZĘ 


were equal to zero, that is, if the center of gravity of each stiffener 
cross section were at the middle surface of the plate, equation (A25) 
would reduce to equation (A3) of reference 9. 


Minimization of equation (A25) with respect to the coefficients an 


yields 
2 2 2 2 
" INS) a (A) - 
|= py DE Mr iba = 
yo\ m xp m 
2- e EL 
EL d o.t/fa/d| A | 
SS - gt a/a) —8 an sin rui. sin nri + 
dD è „p Vm/ dilo — i-i N 
> 2 N-1 
2 EAs%g l Fim sin TH 
N aD 2 d. a N 
i=l 
(EX) RRB, 
F 00 
2 E >, an sin REL = 0 (A26) 
a, | /mx = | n=l 
9 (==) EAgZg 


(EST 2. e. e œ) 
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When each coefficient Fim is determined in terms of the coef- 


ficients a, of the deflection function by the analysis of the pre- 
ceding section and is substituted into equations (A26), the equations 
may then be manipulated to yield stability criterions for the stiffened 
plate. Since the expressions for Fim differ for plates with different 


numbers of stiffeners, separate analyses are made in the following sec- 
tions for a plate with one, two, three, or infinitely many stiffeners. 


At this point it can be seen that the assumption of sinusoidal 
buckling in the longitudinal direction is justified. If sinusoidal 
buckling had not been assumed and an infinite number of longitudinal 
sine terms had been taken in equation (All) for the deflection function 
(and, correspondingly, an infinite number of cosine terms in eqs. (A9) 
for the shear forces induced at each contact line), the terms corre- 
sponding to each sinusoidal component in the longitudinal direction would 
uncouple when the potential-energy expression was minimized and there 
would be obtained a different set of equations (A26) for each value of m. 
This uncoupling occurs only for the particular boundary condition chosen 
in the present paper; that is, points in the transverse edges of the 
plate are free to move in the y-direction before buckling occurs but are 
completely restrained from further movement after buckling. 


Plate with one stiffener.- For a plate with one stiffener there is 


only one shear force induced in the plate - the one along the center 
line: 


F4 = Fim COs = (A27) 


The middle-surface displacements in the plate at the contact line 
due to this shear force are found from the appropriate equations for 
bay 1 and bay 2 derived in appendix B (eqs. (B13) or (B14)), with y 
and k set equal to zero and b replaced by 2d, as 


2 F 
z a lm mrX 
(1) 8 = zo (>) Hat COS e (A28) 
y=0 | 
where 


2 
Zmd 2 ( mxd 2 
(3- u)( + H)cosh == + 2(1 + u) ES + 5 - 99 * u^ Lid 
sinh fusi , mnd a 


a 
(A29) 


1 
a 


ica ———————————À < œ< 
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and the displacements in the stiffener at the contact line are, from 
equations (A10) and (B24), 


00 
= ox) a V Fim MIX , mx < > nx mix 
By” E En EAs B zd H 2 a 
(A30) 


The condition that these two displacements are equal can then be written 
as 


d + Zo] SL = > a, sin = (A31) 
(SE) Mežs ™ 
from which 
> ap sin X 
jm Bel (A32) 


The substitution of equation (A32) into equations (A26), with N 
equal to 2, then yields 


2 
Den 
1 afela) MCA (lay 
1+= — - = + 
2- 2 co 
EI Z d. o.t A 
1 a/a) 7B| a, sin ZE gin BE = 0 
S dD np `” dt | teil 
l + 22] — 


(A33) 


(lo lo ex a 500) 


Equations (A33) are practically identical with equations (Ah) of 
reference 9 for N = 2, when those equations are divided by me, the only 


E ———— er i A cmo IN Rega A 
- ——— A A A E --- 
ee mm asi e a ei m i ae S i i á—À 
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difference being that the flexural-stiffness ratio (called y in ref. 9) 
of equations (Al) is replaced in equations (A33) by an effective- 
flexural-stiffness ratio 


- 2 
Elere _ Eig , 1 EAgzd 
aD aD Ag a 
l + — 
221 3i. 


(A34) 


The criterions for buckling under uniform compression of a simply sup- 
ported plate with one stiffener can therefore be written immediately 
from the results of reference 9 as follows: 


For buckling of the plate with deflection of the stiffener, 


D 
y (ca) d^ | 
= o(m = 2 
EL EAcz m d c.t A 
mu QNIN AA 0 | xD CH o (8/2 a) a) (A35) 
14% = A AA 
ol 9 9 
2 1 
where 


04 = (A36a) 
0o = (A36b) 
and, for buckling of the plate with a node at the 'Stiffener, 
2- 2 
à o.t 
El = |B +p? us (A37) 
£p wa "om 
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where p is the number of half-waves in the transverse direction in 
each bay. The lowest values of the buckling-stress parameter given by 
equation (A37) are obtained when p is set equal to 1. 


Plate with two stiffeners.- For a plate with two stiffeners, two 
Shear forces are induced in the plate: 


Along the contact line at =-$, 
FQ = Fip cos — (A38a) 
and, along the contact line at y= g, 
F, =F, cos REX A38b 
2 "Om " (A38b) 


The middle-surface displacements in the plate at each contact line 
due to these shear forces are obtained as follows: The shear force Fi 


alone produces displacements uj' along the contact line y - -5 which 
may be obtained from the appropriate equation for bay 1 and bay 2 of 


appendix B (eqs. (B13) or (B14+)), with K ==> =-3 2, and b 
replaced by 3d, as 
(u 9 = W +Z i= cog Mix (A39a) 
1 y=-4 Pu) (731 * 432)Rat a 


and displacements up' along the contact line y= g which may be 


obtained from the displacement equation for bay 2 (eq. (B14)), with 


K =i y= 3% and b replaced by 3d, as 


heri m a we um 
J-5 
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where . 


+ cosh sinh —— 


(3 -p)Q + l.) (cosh mad 2d z cane m n 


e 2 (mxd 2 mid _ 
5 - 21 + p^ + 2(1 + pu) (pe ) |cosn a 


2 
(1 + p)? MA sinn AMA 4 5 S +p? AL + p) (m) 


Z3] = l mid 
s sinh 3mxü t mrd 
a a 
(Ahoa) 
and 
3mxd 2mxd _ Zmnd mad 
(3 - wu) + 1) (cosh "Ada cosh x i sinh A )* 
2 
E - 2u + p? + 2(1 + ny cosh md + 
2 t 
(1 + u)? ma. sinh Zad | - 2u + p? + h(1 + ¡ay , 
l dà 
c rc 


nh mrd _ 3mnd 
a a 
(A40b) 


Similarly, the shear force Fo produces displacements uj," along the 
contact line y =-5 


" _ Ll 
(uy pint = ES (231 x ZA: cos == (Alla) 
and displacements up' along the contact line y= 5 
( 1j = yay (Z + z ) em cog MX (A41b) R 
"2 -4 2\mn 3l 30/Edt a 


"2 | 
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P4 


Ihe total displacements u; along each contact line are obtained by 
adding the displacements u;' and wu," due to each of the forces and 
are 


e F F 
(m$ = 5 2) (ts + E30) Ra (231 - Saal -i COS a (Ah2a) 


(A43b) 


The conditions that the contact-line displacements in the plate and 
ín the stiffeners are equal can be written as 


A F = 
- +3 (23) - 242)-8 — £n... 5 a, sin na 


1+ Z (Z + 2 
2 dt ESTA ESE n=1 


o. EA 


(Allia) 
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F F = > 
=(Z |= Z JE z > noe. + | L+ 2(Z3] + 23) 3 —— iż = an sin ani | 
21/31" “32/at „3, — 2 D ee 3 
ru) EAgZg (Ex) EAg2g | 
( Alib) 
from which the following relations are obtained: 
F co | 
im == l ap (sin 77 + sin ZE) + 
std z zl 
(Jasa 1 + Z5 = n 
L l A; » a, (sin DX _ gin a) (A45a) 
=1 
1 + 230 at | 
F co | 
em 1... 1 im uM 
3 aż 2_ %n(sin E + sin ES) | 


Substitution of equations (A45) into equations (A26), with N 
equal to 3, then yields the following groups of equations: 


For buckling symmetrical about the plate center line and with deflection 
of the stiffeners, 


at Es (*68+1 p 2.6545) s 
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= 2 La 2 00 
I BAaż d o.,t Ą 
xS + 1 — Teis - Se la) - Pm ET )-0 
Dos, Ag @D xp \m/ dt] ¿y | O8+1 B+) 
a at 
(A46b) 
(P=0, 1, 9, 3.440) 


for buckling antisymmetrical about the plate center line and with deflec- 
tion of the stiffeners, 


Lon PP NN: EAN 1 eee 


A a 2 8 
i EAgZzą d ot a/d Ag B 
E prt Rn (248) at 2. esa? sgh sO 


s=0 


- 2 C= é 
EAgZzg AB 2 E " a "" 
As aD ER M dtl eso ( 6842 6844) 
dt 


(A47b) 
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and, for buckling with a node at each stiffener, 


2 
TEM E ba = 0 (448) 


xD 


Equations (Ah6) to (A48) are similar to those into which equations (Ah) of 
reference 9 for N = 3 separate. The differences between the equations 
of the present paper and those of reference 9 are that, instead of the 
flexural-stiffness ratio 7, equations (A46) contain an effective- 
flexural-stiffness ratio 


o: BER (Ah9) 


and equations (Ak7) contain an effective-flexural-stiffness ratio 


S 


- 2 
eff _ Elg 2 1 EAgZg 


dD 
1 + Zza 3 


(A50) 


E 


The stability criterions can be written immediately from the results of 
reference 9 as follows: 


For buckling symmetrical about the plate center line with deflection of 
the stiffeners, | 


2 22 
a eee EAgZg a RAM xD E e Se “s 
Ag aD sin 04 sinh 05 xD m dt 
L L Z — 
31 dt 01 05 
V3 V3 
a cos 01 zw cosh 0 


(A51) 
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for buckling antisymmetrical about the plate center line with deflection 
of the stiffeners, 


2p sala) at " o 
—- l EÀsZS _ zm | xD yd Gxt (a/a) Ag 


dD sinh 05 sin 07 £D m dt 
l + 230 at 785 — E 
dt | o E 04 
JE + cosh 05 G + cos 04 
(A52) 
and, for buckling with a node at each stiffener, 
Da 
d^o, t e 
x = p? 2/0 (A53) 
xp la/a m 


(P=1, 2, . . . 0) 


Plate with three stiffeners.- For a plate with three stiffeners ? 
three shear forces must be considered: 


Along the contact line y = -d, 


F, = F, cos = (A5la) 


Fo = Fy, cos — (A5Hb) 
and, along the contact line y - d, 
F3 = Fz, cos SIX (A54c) 


The middle-surface displacements in the plate at each contact line 
are obtained from the displacement equations for bay 1 and bay 2 derived 
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in appendix B (eqs. (B13) and (B14)) by the procedure indicated in the 
preceding section as 


e | Fin Es 
io E (=) Bos. 3 Ml + 3m 232 T +0 lm ~ Ao) la z mX 


(A55a) 


= (E) |> ha, Be Eat T Mi E + 3m z os um a” (A55b) 


Fon , i Fan 
DAŁ as = ax) (20) Bess ~ on) t + "moe + 5n + hom) cos 2 


(A55c) 
where 
(3 -p)(1 + u) (cosh “ant Amd + cosh ma. z anne sinh aa) + 
2 
b- OL + u? + 2(1 + (est oa EDI, mM 
2 
2(1 + u)” Pd sinh Sh + 5 - S + a? + 6(1 + p) (maa) 
Laa = 1 md 
dis sinh imę + ar 
(A56a) 
(3 - H)(1 + u)(cosh "md - cosh Emad E Eur sinh Ems) + 
2 
5- 9n « pe + 2(1 + (ust | cost end + 
| 2 
2(1 + u) zx sinh md z - +p? 4 6(l + (a? 
_ i md | 
Aom ^ a 


ah uen _ lmnd 


a. 


(A56b) 


a wać - — 
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(3 = o + w) (com 3 IA simn m) «| 5 = a + us 


2 
h(1 + y co a mid _ (1 + u) mae sinh um 


Mam = p OE 
à a sinh rd, Mid 
(A56c) 
Imra 2/mxa Ye 2 
1 mą (27 HG + u)cosh "ES + B(1 + u) (Tae ) +5-%+u 
Mun === ———— 
Y a WIR 


(A56à) 


The displacements in the stiffeners are 


2 F co 
z -( a lm mix , mt > > nn micz 
) EAS COBRE 28 - a, sin P cos "= 


z (Sx Ju - (2) al i Ea, sin = cos = 


(A57c) 


53 


-mr m ee —- = m —À — — ee —- 0 
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The conditions that the contact-line displacements in the plate and in 
the stiffeners are equal then can be written as 


A F 2 
(sim - Lal S — 2 > sin T (A58a) 


A F A F 
dom qe = + (1 + Ayn) = z* 
(E) EAsżs BAgig (=) 
Ag F = nx 
tag, E (A58b) 
=) FS 
1 As ^ Fim Fon 
504m 8 hom) ae cr Xam - Gum 
Cz) zg (=) EAgZg 
1+ (Mm + on) oe ELEC SN = x a, sin 2% (A58c) 
nt TW n=l . 3 
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from which are obtained 


Fim oy gE la E sia 2). 
Tu, Ge R] + l-a gs) T 
“3m dt 5 > an sin PL + 
=1 
(2 x EG + Z): sha E) 
no 2 lnn E - sin 322) (A592) 
dt 
“2m == a Z >> a (sin 22 + ain 32%) 4 
(25) "Engi (a + om 2) (a + Dym B) - SU =} = 
l + Am PESE. (A59b) 
(2 + ya ZI * Mm G8) - sha Z) 
m --i|- L + Mm zę 2a (sn BE + sin 320) , 
(STs, (a + in 8) Q + Ma = E SU ss) iin 
Ds at > en sin BE 
(2 + dam a). + Mm ac) - Da ac) 
EE >. a, (sin nx - sin 3r) | (A59c) 
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After equations -(A59) are substituted into equations (A26) the 


following groups of equations are obtained: 


For buckling symmetrical about the longitudinal center line of the plate, 


Elg | 1 + (Pm + Mm ym, Bs SAN. 
aD 
(2 + Man ZR + m GE) - a 38) 
2- 2 co 
Ea) E 2 (28841 ~ 89.547) ü 
Mm > Mim Ag - 2 
E db EAgZg 


y dD 2 (^843 ~ *85+5) = ? 


afs t/a/a\? Ag| © 
= (2/2) dt = (38841 9 8547) N 
(Ma 5 "im ) As _2 
SE 2 pias 


2 aD UN (?8243 7 88545) 7 U 


Lr =0, 1l, 2,.. 


(A60a) 


. œ) 


(A60b) 
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2= 2 
a + ( + [2r 4 > (2/2 ka `. Sixt (a) a + 
El 2 Y 8r+3 


Je lo Me + 84 a 
EA A et 
ZIEL 


a°G,t [a /a V? Ag | Y 
(e) : E (^8s«3 7 *8s45) 7 


Cua; eye 


-2 c 
dt EAgZg 
o aD 2. (88+1 m 2954-7) 59 (A60c) 
Ag Ag Ag is 
d + Alm ON + Mm E) a sla =) 


(r=0, 1,2, ... œ) 


: L + (m ds Vinis EAgZg? 
2 aD 
(2 + dam ZIR + m 32) - sha Z 


a“ozt (a/aYé Ag| Y 
EE) : 2 (98543 ği 28845) z 


| (un = Men ki ia) As | r" 
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for buckling antisymmetrical about the longitudinal center line of the 
plate, 


t 2 
Elg | 1 EAgźg” — d'Gyt (2/2) 
Ag dD xD m 


A co 
at UM ( 0842 6346) 
em dt 


b+ er FEST- MTL. 


E, 01 ETA 


at | o (8542 7 98546) = 0 


and, for buckling with a node along each stiffener, 


poto e 
[5 2 1. 2, . . . œ) 


A comparison of equations (A60) anà (A61) with equations (Ah) of 
reference 9 indicates that equations (A61) differ from the corresponding 
equations (Ah) only by the substitution of an effective-flexural- 
stiffness ratio 


-—————— M ——— AM ————— — —— — z — - 
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(A63) 


for the flexural-stiffness ratio y. Equations (A60), however, differ 
radically from the corresponding equations (Ah) of reference 9 in that 
those equations containing deflection-function coefficients of the form 
888+1 and. 2954-7 are not independent of those containing deflection- 


function coefficients of the form a and a - A single criterion 
88+3 88+5 


for buckling symmetrical about the longitudinal center line of the plate 
is thus obtained rather than two criterions corresponding to those of 
reference 9. Application of the method of solution of reference 11 and 
the results of reference 9 yields the following criterions: 

P4 


For buckling symmetrical about the longitudinal center line of the plate, 


3 |a2> 
Ñ p = in 8 inh @ 
as E 
| < " 


rod cos 61 rn cosh 85 


Mua me, ad maż 3 MTS, 
T sinh 9 sin 8 
TES E 
Va v2 


= + cosh 95 EJ + CO8 04 


=0 (A64) 
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for buckling antisymmetrical about the longitudinal center line of the 
plate, 


» Alla? fe ext 


- m 2- 2 
CMT NE. RA E O Ag (165) 
* Mm gt 85 F 


and, for buckling with a longitudinal node at each stiffener, 


2 2 
d Gyt _ o a/a s M 
ŻE (+ s/a) G OM MET. (A66) 


It is noted that, when Zg is equal to zero, equation (A64) yields 
the two stability criterions for symmetrical buckling given in refer- 


ence 9, that is, equations (AT) of that reference for = = and 2 = 


= 
U N 


ElL 


Plate with infinitely many stiffeners.- As the number of stiffeners 
is increased the work involved in obtaining stability criterions becomes 


greater and the stability criterions become very much more complex. A 
complete investigation of all the possible modes of buckling of a plate 
with infinitely many stiffeners is therefore out of the question. The 
investigations in this section are limited to that mode of buckling in 
which the deflected shape of each bay of the plate is identical with 
that of the others, since the numerical results of reference 9 indicate 
that this mode of buckling is predominant for a plate with infinitely 
many stiffeners. The deflected shape of the plate in the transverse 
direction is considered to be of the form shown in figure 5(a), sym- 
metrical about the center line of each bay and horizontal tangents at 
each stiffener. The equation of the deflection surface can then be 
expressed as 


m 
w = sin T > an COS E (A67) 
n=0 


The shear forces induced after buckling are identical at each contact 
line and are distributed between adjacent bays of the plate as shown in 
figure 5(b), half of the force being applied to one bay and the other 
half to the adjacent bay. 
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Only part of the buckled plate, the part between the center line 
of one bay and the center line of an adjacent bay, need be considered 
since the deflected shape is periodic. The potential energy of this 
part of the plate can be written as 


à 2 2 
9 Rr (v2) «21-04 dw 


O 5 òx òy Ox dy 
E 2 
=, B a dica JN > (ów 
SEE ABE = Fez == dx = 
aD Jo Max o ^b d - 
z - 8 z 2 
Gxt [ 2 (Qu Bls J “(Ow 
= Jo MC) odc 0 t3 = 282) 


where EK, is the shear force induced after buckling. 


If E is taken in the form 


FP .=F! cos m (A69) 


the displacements u,' in the plate at the contact line are derived in 
appendix B (see eq. (B20)) as 


vey (3-u)simh REL (+ wR EN — 
Uu) 773 à mt gc COP a (ATO) 
=0 cosh HE - 1 


— mr mr Luca ee AA ——M 
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and the displacements Ug! in the stiffener at the contact line are 


= "EX zoe 
Ug! + (as) - -(&) TW cos ral x Za 2 pa cos MÆ (471) 


The condition that these two displacements are equal yields 


Ft Sw 
eee AT2 
- (A72) 


und 
-— _ (3 - p)sinh —— - (1 + ij A 
(m) 1,131 a m48 


4 PL 1 a dt 


Substitution of equations (A67) and (A72) into the potential-energy 
expression (eq. (A68)) and minimization of the resultant expression with 
respect to the deflection-function coefficients a then yields 


the equations 


2 


(e) G + oro) « Se’ 


aD mid. aD — 
1 + y 3 = u)sinh —— - (1 + TR == a mad, ^g 


P" agi » mE a dt 


d o,tfa/a| A co 
= e/a) E w (r=0, 1, 2, . . . c) (A73) 
nz 


where 5,9 is the Kronecker delta (Bro =1 if r=0 and Ba =0 
if ro). 
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When the method of reference 11 is used to manipulate equa- 
tions (A73) and the results of reference 9 are used, the following 
stability criterion is obtained: 


- 2 
dD mrd mad dD 
- sinh — ~- (1 + pacta 
DETAIL u) 7 (1+H)— mua AS 
4 cosh =: 1 E 
a^, t a/a 2 Ag 


This equation is, except for the flexural-stiffness ratio y being 
replaced by an effective-flexural-stiffness ratio 


2 
dD dD mud md dD 
(„lu (3 - k)sinh “= - (1 + M7 mra Ag 
h cosh TAL _ 1 a di 
a 
(A75) 


identical to equation (A7) of reference 9 for a equal to zero. 


The terms in the deflection function that yield the stability 
criterion for buckling with longitudinal nodes along each stiffener 
have not been included in the investigation of plates with infinitely 
many stiffeners. This criterion should be considered, however, and is 


given as before by 
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APPENDIX B 


MIDDLE-SURFACE FORCES AND DISPLACEMENTS DUE TO SHEAR FORCES 
APPLIED ALONG A LONGITUDINAL LINE IN THE PLATE 


Plate of Finite Width 


Consider a plate subjected to an arbitrarily placed middle-surface 
shearing force of the form 


F = F, cos — (BL) 


along the line y=k 2- (See fig. 6.) This force divides the plate 


into two bays. Within each bay of the plate, middle-surface forces and 
displacements are continuous and are related by equations (Al) to (A3). 
At the line of application of the shear force, however, middle-surface 
shear forces are discontinuous and forces and displacements of the two 
bays are related by equations (A5). It is convenient to introduce a 
force function Ø; for each bay, such that 


2 
dØ 
Ny = jm A 
Jm dy 
2 
à 0 
Ny L (B2) 
à 
i (3 = 1, 2) 
37% an 
XY sm 0x dy 
and 
vig, = 0 (3 =1, 2) (B3) 


p e A —— —— 
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These force functions satisfy the equations of equilibrium in each bay 
(eqs. (A1)) and ensure the compatibility of the strains derived from 
the middle-surface forces. Suitable stress functions for the present 
problem are j 


NE (An + Bym) sima 7 + (Cin + = Djn)cosh sin = (8h) 


From equations (B2) and the stress-strain relations given by equa- 
tions (A2) and (A3), the following expressions for the forces and dis- 
placements in each bay can be obtained: 


For the forces, 


2 r 
mty mxy | ¿mx MICX 
(Ca t Bam + EC D a )cosi EUN (=) sin E (B5a) 
- nu my 
len (4n + a Bjg,) sinh a + 
muy Inxy |/ mna 3 micz 
muy 
NY jm - (An + D im + Eu Bjm) cosh E t 


i 2 
(C + Bam + — D jm) sinh aa (=) cos =x (B5c) 
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and, for the displacements, 


Etujm = - la + à) (An $ = Bim) + 2 simm ay + 


|o + H) (Cm + =v Dj) + 2 cosh ay = cos ma (Ba) 


EtV im m oce |a + u) (Am + — Bim) - (1 - 1)D yu cook = -+ 


[a + w(eja + 2% Dyn) = ©- vR Q sin BS 


(B6b) 


Rigid body displacements are left out of the expressions for Uim and 
Y jm since they do not influence the solution of the elasticity problem. 


It is seen that the forces and displacements already satisfy some 
of the boundary conditions that are imposed. The displacements in the 
y-direction and the normal forces vanish at the edges (x = 0 and x= a). 


The conditions at the line of application (y = Kb) of the shear force 


(eqs. (A5)) yield the following conditions that must be satisfied by the 
unknown coefficients of the stress functions: 


(^om - A44, )cosh D t (Bom - m Cc cosh D + Binh m + 


kmstb kmnb kmab kmnb| ./& 
(Con - Gia) 8inh a + (Don - ndl Da sinh zm + cosh uo) | ) Fm 
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. kmrb kmsb kmnb 
(Ap, — Ayn) sinh "uli (Bon - Bin) > inh —— + 
Km D kmzb kmrb _ 
(Con - Cy 7) 95H + (Don - Dim 2. osh = 0 (B7b) 
kmsb kmnb kmstb 2 b 
(Ao, — A,,,)8inb da + (Bon - Br) sinh Toa + Tou cosh j + 
kmrb kmxb kmanb 2 kmrbY _ 
(Con - Cy )cosh Da + (Dom = nll zm cosh m T gem sinh zm = 
(BTc) 
kmxb l kmaub kmtb 1-H kt | 
(^on Az q) cosh oe (Bom Psy a.  COBh "go = RH sinh ==) + 
Km D kmnb kmtb 1-1 kmxb\ _ 
(Com ~ C),,)8inh Oa + (Dom - ndl Sa sinh "Ba - lud cosh Da ) = 0 
(RTA) 
From equations (B7) the following equations can be obtained: 
i o l+pfl-u kmxb . kmrb kmsb |/ a 
A edP T 6 fi cosh A ZA sinh ZW X: ) Fn (B8a) 
1 + kmstb 
Bom - Bim = æ x ginh EE ==) m (B8b) 
 13+p4/f1-p4 kmnb , kmxb kmstb\ a 
Com - Cim > G ng m cosh A X ) Fin (Boc) 
1 + kmstb/ a 
Dom - Dg = = cosh S (=) m (BB) 


- o. -- MÀ Án — — X om 
a aa See ee ip ee cg my aa m A ATT nn a, a oig — ~~ 
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Four more equations from which the coefficients of the force functions 
can be determined are obtained from the conditions at the edges y = 2 


of the plate; that is, the middle-surface forces normal and tangential 
to these edges vanish. Then 


A, Sinh ETÈ + p, Mtb sinh ED + C cosh DAD + Dom B=? cosh BED = 0 


em 2a em 2a 2a 2m 
(B9a) 
A, cosh mtb + Bon (Be cosh Ja + sinh zab) " 
: mrb mrb mrb mrb 
Co, sinh mde: nal sinh > cosh > ) = 0 (B9b) 


A, sinh B® - pg MD ginh WE | C. cosh BM + p mb cosh B® -= 0 
lm 2a lm 2a 2a lm 2a lm 2a 2a 
(B9c) 
. A. cosh WD _ p (MD cosh WÈ + sinh Mb) _ 
im 2a In(2a 2a Da 
b mstb mb 
C... sinh ED + p (Z2 sinn E? + cosh BP) = 0 Bod 
im 2a 1m| 2a Da Da (B94) 


from which are obtained 


-a(g sinh ox + 


kmstb mnb mxb kmxb 
% EDI cosh(2 - k) + 2 cosh dB 
L = U mstb l mb kmstb 
= |sinh(2 + k)—— - sinh(2 - xk) + 2 sinh = 
: > Le A os po) 
Aom + Am 4 (5) Em mtb  mxb L ) 
sinh ——= - —— 


ee G. a 
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cosh(2 + En + cosh(2 - DS + 


3-1 kmb  bkkmaxb imb 
: 2 nm cosh Da + Da sinh Em 
l L U 2) 
Bo, + B = = —__ fy — F. 
lm = m ; ~ 
4 sinh WP mxb 
a a 
(BLOb) 
l- pu 


mab mrb 
cosh(2 + k)—— + cosh(2 - k)— — 
Plese + x (2 - wje 


kmstb kmstb mrb 
2 cosh + sinh(2 + k)—— - 
2a ea p ( IS 


sinh(2 - k)4ID _ 2 sinn Emxb| . 
( JR zm 


2 
y (=) GS kmnb 


1 2 ca 2a 
L T Ha 
C +C = —) PF 
em — "lm h (riz) i sinh EXb , mb 
a a 
(BLOc) 
mxb MTD 
sinh(2 + k)5 - sinh(2 - k)5 2 - 
2 3 -H sinn Emxb _ Hmb sogh Eb 
2 lt ca 2a 2a 
Do + H - ithe) P 
em lm K Jna a giną WÈ _ mtb 


(BLOQ) 


a — — — — — 
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Solution of equations (B8) and (BLO) then yields 


2a, ca 


2a 


(2 l-g, an cosh SB _ lac = x (mee E 


l- HE lai nn Emxb , Eb cosh(2 - k W 
2a 2a 2a 


l +y 
l -p mstb 
2 — k)—— 
2 l + sinh( p 
A. ae (zz) F 
lm m. 
n sinh IP _ mrb 
a a 
(Blla) 
-p kmxb 
cosh - 
L +L U 2a 
o 2(i- k)EXP sinn ED + cosh(2 - k) 
ca ea 2a 
Bim = = ux) Em 
x 
4 sinh mrb " mnb 
a a 
(B11b) 


o LJ 
mrb l - kmstb 
lec > x) (52) z LI kom x 


(2 Ica S sinh SÓD + 
L +L U 2a, 2a 


l-u mxb kmxb mxb 
h(2 - k)—— - inh(2 - k) 
- o > Cosh a żę sinh( a 
H/a 
Cin =-— (ix) Em sinh Mb, mb 


(Bllc) 
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kmab 


mrb 
2(1 - k)— cosh 
( a cos 


l J - M ginh Em"b _ ginh(2 - «Bt 
L +L U 2a 2a 


Din = 2 (2) R —————————————————-—— .— (mua) 


and 


2 
l - pAimxb kmxb mrb 
(x En a cosh zs E t E ) 


I- E san kmrb + TT cosh(2 + kj + 


1+pq4 2a 
l- yu : mxb 
2 > 
: 5 BEST sinh(2 + KS 
— + H/2. 
Am h (mr) Em sinh wtb _ mb 
a a 
(Bl2a) 
j 7 V cogn E 4 2(1 + Hb STB 4 
1 L U 2a 2a 
cosh(2 + gn 
l + pya 
Ba =- —) E 
2m R m 
* sinh ma + => 
(B12b) 
mrb : L - H | kmnb 
2(1 + x) (Ba) ane cosh 9A + 
l- u mrb ( l-u 
cosh(2 + k)—— + [2 = 
1+ ( E» 1+u 
x) mab sinh kmnb + kmsb sinh(? " k) mtb . 
l + usa 2a ca ca 2a 
Com das m (ax (B12c) 


A a OR PE m A ai IN TT A —— M M AM me 
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mrb kmrb 3-H . kmxb 
Dom = l El En mrb mrb REM 
Bini RN 


Of particular interest are the displacements in the x-direction of 
the plate middle surface. These displacements are determined for bay l 
by substituting equations (B11) into equation (Bóa) to obtain 


= (3 - n) +p) cosh(2 i "JE vS + cosh(2 _ ey | des x 


um "| sinh BAD y zb > 


| - 21 + pe + 2(1 + y) “(1 - v, - 2) (22 Loa (e + x) a + 


ERZE OCZ 


3-00 + A E yz |. 


sinh z a eae 
a 


a 
5-2 + 2(1 +n) - x) - 21) een 5 + 1) Re - 
st ata a Dl col aR 
(1 + u)? ze (a - w)stnh(2, + = a = + (8 + x) sina (2 = ey = SE E 


(3 - H)(1 + age mp x) sinn (= + x) ae + 


(2 + - vem - Jem TT cos mX (B13) 


o 
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and for bay 2 by substituting equations (B12) into equation (B6a) to 
Obtain 


= a A = 2y mb ey mb 
U2m = mxb . zł H)(1 + 1) [som BE NE + cosh(2 Pus + ga + 
a a ` 


b - Mapes NA + u)?n * k) (2 + SETS + yan + 
patata E (er. JE 

(i + p)? snl + x) sim(2 ++ x) Ba a (= z k) sinh(2 š 2 + A + 

(3 - 4)(1 + yn (2 + F + x) sina (Er + x) a - (2 - X + vlan E zi 5 


1 2y )zn ( 2y I 
— - 1 + osh| 2 — + k|— — 2a c —— | = 
„akg | 1) ( 1) cosh( + S = cosh = + kJz- 

a a 


b- OL + y? + 2(1 + u)?(1 * as * eros IZE + x) Bm + 
b- a +u? + 2(1 + u)*Q + KI, - E Tag - VBR 
(1 + p)? B lS, daal + ZE (22 doal - + ae 


(3 -p)(1 + We (2 + RU + i) sinn (SY + x) + 


(2 = zi + x) sinn (= - e] 5 = = cos = = 


AH mc M ne RM Z r Piba PP D w Pe m na 
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Plate With Infinitely Many Stiffeners 


For the assumed buckling configuration of a plate with infinitely 
many stiffeners, a separate derivation of the middle-surface forces and 
displacements must be made. It is apparent from considerations of 
symmetry that only one-half of a bay of the plate need be considered and 
that the conditions that must be satisfied at the edges of this half-bay 
are given by 


zy! = = 2 
Ny =v O at y 5 
W a at y=0 | (B15) 
Xy 2 
YT =0 at y=0 
The shear force F, may be expressed as 
Fo = Fp’ cos EZ (B16) 


and the appropriate force function for the half-bay is given by 


t- t WY as ated t DY pt EU > 
ba ac de B')sinh > + (Cy > Dy”) cosh = sin > (B17) 


which is similar to equation (Bh). Expressions similar to equations (B5) 
and (B6) are obtained for middle-surface forces and displacements in the 
plate. The substitution of these expressions into equations (B15) yields 


| (1 + pap! - (1 - H)Dy' =0 (8188) 


A! + Dg! = (=) R. (B18b) 


mx 
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A y Pd o, (z mid L < U S 
hy cosh 7 By (5 cosh Z low e + 


MIA. mna. mn. l- pu mid 
Ç inh — + Dy (ES sinn EE - AS) = 0 B18 
Cp" 8 2a 2a 2a Ju te 2a (B18c) 


AAT cosh S + B! e cosh zu + sinh y + 


md mrd mad. mxd \ _ 
Cn’ sinh Da + Dy'( 3 sinh Da + cosh mad | = 0 (B184) 


The solution of these equations for An's B,', Ca: and D,* yields 


2 

ia 

A = HË) Fut (B19a) 
ME goth mada Y P | (B19b) 
E: 4 2a Mx m 
l-H ginh Ad _ mad 

L L N L +L U a a ją = 
ARE E rcc AA C F * 
œ 4 cosh P _ 1 (Gz) > (819c) 

a 

AE (B19d) 
m "y unt m od. 


The displacements along the line of application of the force, which 
are of particular interest, are then given by 


o Ge Main FE TREE, R 
(Um - ux wol C GACER O (B20) 
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Forces and Displacements in the Stiffeners 


To each stiffener, forces per unit length of the form 
Fi; = Fy, cos sa (11, 2, v (WS KR (B21) 


are applied. From equations (A6), the force in the stiffener is related 
to the applied shear force by 


IPL = p cos MEX (1=1,2, . . . (N- 1)) (B22) 
a 


in which case 


Py = > Fip Sin AT (1=1,2, . . . (N- 1)) (B23) 


Constants of integration vanish since each load P, must vanish at the 


ends of'the stiffeners. Displacements in the stiffeners are then 
obtained from the forces by the use of equations (A8) and (A9), which 
yield 


e F 
ug, = (E) BA, cos EE (1=1,2, . . . (N- 1)) (B24) 


Constants of integration are omitted since they represent rigid body 
motions and have no effect on the problen. 


+ A 
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TABLE 1 


VARIATION OF Zyq WITH BUCKLE ASPECT RATIO 


L. 
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TABLE II 
VARIATION WITH BUCKLE ASPECT RATIO OF mn FUNCTIONS FOR 
BUCKLING SYMMETRICAL ABOUT THE LONGITUDINAL CENTER 


LINE OF A PLATE WITH THREE STIFFENERS 


SHO 
MGMT 


ON O co 
QD ed an 
CONO “A 


1. 
1. 
1. 
Le 
1. 
1. 644 
1. 782 
1. 
2. 
2. 
2. 


29 
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Figure 1.- Compressed plate with equal and equally spaced longitudinal 


stiffeners on one side. 
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Figure 2.- Hypothetical extension of stiffeners to contact line. 


— A am —— SO ie —— A — — — 2 


e | ——— — -— 


Figure 3.- Functions appearing in expression for effective flexural 
stiffness of stiffeners. 


CAg2 NL VOVN 


(a) n1» Mar Age and ls” 


Figure 4.- Functions appearing in criterion for symmetrical buckling of 
plate with three stiffeners. 


£lgo Mb VOVN 


E9 


eaid -— x 


(b) n, 8nd mg. 


Figure 4.- Concluded. 


€)9¢ NL VOVN 
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(a) Shape of deflection surface in y-direction. 


(b) Shear forces at contact line. 


Figure 5.- Conditions assumed for plate with infinitely many stiffeners. 
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Bay 2 


TT. US 


F= E cos M7 X 


Bay | 


psu "Mese 


Figure 6.- Shear force applied &long & longitudinal line within the plste. 
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